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14 . Aktiraaf * 

Numerical and experimental results are presented on the pullout 
phenomenon in composite materials at a high rate of loading. 

The finite element method was used, taking into account the 
existence of a virtual shear deformation layer as the interface 
between fiber and matrix. Experimental results agree well with 
those obtained by the finite element method. Numerical results 
show that the interlaminar shear stress is time-dependent, in 
addition, it is shown to depend on the applied load-time his- 
tory. Under step pulse loading, the interlaminar shear stress 
fluctuates, finally decaying to its value under static loading. 
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1. Infcr'oduetlon 
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The pull-out phenomenon of flbui's from a matrix 1b an 
important problem closely related to the stress transmission 
mochanlam and fiber stripping processes in composite materials. 
This phenomenon is also a crucial factor of the impulsively 
absorbed energy, for example, obtained from impulse-teats on 
isots tljRJ, as has been already pointed ouv- by Outwater [33. 
Lawren'^o investigated [*]] the pull-out phenomenon of fibers from 
matrix Iviyers under static loading. However, for a high rate of 
loading, it is often insufficient to rely only on the usual anal- 
ysis of static loading. Namely, in static analysis shear stress 
waves propagating within a material have not been considered. 
Therefore, we investigated the pull-out phenomenon of fibers from 
•'omposite materials at high rates of loading, numerically and 
experimentally, taking into account propagation of stress waves 
within materials and existence of a virtual shear deformation 
layer as an interface between a fiber and the matrix. 

S . A pull-out model and equations of motion 

Let us consider a one-dlnensional model of Figure la to 
examine the pull-out phenomenon of fibers from fiber-enforced 
composite materials (or piled boards). In tbe model, two differ- 
ent layers are connected by a massless shear deformation layer 
(an adhesive layer). Through this shear layer propagates the 
shear stress that is proportional to relative displacement and 
inversely pi^oportlonal to the layer thickness. The shear layer 
(the adhesive layer) in the model can be a matrix layer, or a 
kind of shear stress px-opagatlon layer (transition layer) at 
the interface between fibers and the mati’ix, if one considers 
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Material 2 aa a matrix. In t4io one-dlmenalonal models Figure la /62 
and lb are equivalent to eaoh other, 


Figure 1. Idealised 

pull-out model 



If a fiber la pulled out Impulaively at one end, the load 
propagates through the fiber as a stress wave and reaches a dls- 
eontlnuous junction. For static loading, one can analyse forces 
acting on a fiber-matrix interface mid stress variations within 


the fiber only by examining balance of forces without consider- 
ing the inertia of a material. However, for dynvamteal loading, 
where stress Wvives propagate, one has to consider balance of 
forces taking into account the inertia of the material, that is, 
equations of motion. For simplicity, let us assume that layers 
are flat and two-dimensionality can be ignored and consider pro- 
pagation of stress waves in one dimension. Prom Figure 2 wo see 
that, as a pull-out stress wove reaches a discontinuous junction, 
equations of motion are written as 
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Figure 2. Oomponents 
of forces 


I] V . » S, » ^ 




I*' >T 

* f 2 Ik V 


f. 4 


I*. 

tC*, k k ts jr ft 

"S t sv i 

o; 


Where u, E and p are the displacement, Young’s modulo and the 
density for each material, respectively while the suffix 1 or 

(']) (p) 

li corresponds to the material 1 or 2. il.lso A and are 
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cross sections of materials, given respectively by 
G stands for the stiffness (shear elasticity) of a shear layer. 


3* Numerical analysis 

It is extremely difficult to solve Equations (1) and (a) 
analytically. They have been solved only under limited boundary 
conditions C5] . It seems almost impossible to obtain an analyti- 
cal solution for the pull-out phenomenon, since the boundary 
condition at discontinuity is complex. Therefore, in oi»uer to 
obtain a realistic solution we analysed equations numerically. 
Although any of the differential methods, the method of charac- 
teristic curves or the finite element method based on the 
Hamilton principle, may be applicable to this problem, we have 
employed the finite element method. 

According to the Hamilton's principle. Equations (1) and (2) 
are equivalent to a functional Hamiltonian given by 
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As shown in Figure 3, the triple layer junction of fiber- 
shear layer-matrix is meshed into finite number of elements. Let 
us denote two nodal points of the i-th element as i and i+1 and 
write its length as L. Displacements of the point i at time t 
for both materials 1 and 2 are written as 

if. n. ' ■' , ■*' /'ij ^ 
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Figure 3. Finite 

element mesh 
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For displacements of the point i+1, one writes 
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Assuming linear variations of displacements within an element 
and interpolating by a linear equation of the two nodal points’ 
displacements, one may obtain dlsplacomentiB within the 1-th 
element as 
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Substitute Equations (S) and (7) into Equation (3). Take var- 
iations and set ..// o . By adding an external force on the nodal 
points, one gets the equation of motion for the element 

/■’!«: p. (8) 

where i-q and f/j are the nodal point displacement vector and 
an external force vector at the point, respectively. Their com- 
panies are given by 

{»ts* ft,'*'*, .1,,, . (9) 
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An element mass matrix [ml and the stiffness matrix [kl are 
given respectively by 


,1 >/.*'• » i y'» n 

i> :M' v .1 


>1 ' ! i* 

t ' 

I >1 


( 11 ) 


'*J j 






1 5!-.#^ j .l.iyti- 

I Ui.f 

I a-.. ^ .V 

!• ' * i> 

3 i! , 7 . 

l! / 


Ht, 


r «f.T. j 1 
. , * 


I ili.f J‘S‘. 


I i! «.f 


1 111# I is.j 

It ft .5 ‘ * / 


1 » I.#, 4‘-‘T‘ 


1 III./ 


763 


1 5f« 

:i 


1'.. f";*' 




These are called a eonslobenfe jmss matrix am! a eonsisten., 
stiffness matrix. Instead one may use lumped matrices for siw- 


pllelty. They are given by 
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^1 . Experiments 


We have experimented with tlie impulsive pull-out phenomenon 
using pull-out test specimens . Results will be later compared 
with those from numerical analysis by the finite element method. 

^1.1. Test specimens 


Figure shows the shape and the sine of a teat specimen as 
well as locations of strain gauges attached to it. A shear 
layer is formed by epoxy-adhesives. To Investigate the effect 
due to the thickness of a shear layer, we used two kinds of test 
specimens, one with a 0.05 mm thick shear layer and another with 
1 mm thickness. Both materials 1 and 2 are made of aluminum. 

We used commercially available epoxy-glues (Araldite) as adhesives. 
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Figure E|, ipest specimens 
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Figure 5 . Experimental 
apparatus 


k 2 . Experimental apparatus 

Figure 5 shows the apparatus of our experiments. Let a 
weight free-fall from the height of 700 mm and collide with a 
stopper. The collision generates stress in the rod until the 
weight leaves the stopper. The stress then propagates as a 
stress wave towards a test specimen check. When the stress vmve 
reaches the check, a part of it is transmitted into the Material 
1, while the rest propagates back inside the rod as a reflected 
stress wave. The transmitted wave inside the Material 1 later 
arrives at the adhesive junction with Material 2, and the phen- 
omenon of an impulsive pull-out takes place. Though the shear 
stress at the interface between a fiuer and the matrix is a 
crucial factor in the phenomenon of pull-out of fibers in com- 
pound materials, it is extremely difficult to measure this shear 
stress or shear strain. Instead, we have measured for simplicity 
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bimo variafclons of longitudinal strains fox* Maberlala 1 and S 
using aeeurate strain gauges and compared with numerical results 
gnala from sti^aln gauges were amplified by a DG-amplifler 
(Shinko DA-4007PJ Frequency quality 50kH!3/*3 db) and then 
recorded in digital memories. Strain gauges were S mm long 
(Kyowa KPG-2-11) . 













Figure 6, Stress wave records of A specimen 
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Figure 7. Stress wave records of B specimen 


5 . RgauXfcs and diaousBlona 


5.1. Expenlmentol noaulta 


Figures 6 and 7 show t4tne variation of strains at eaeh 
gauge loeatlon in test apeaimeno A and B measured in the impul- 


sive pull-out test. Here time is set at t » 0 when the front 
of a stress wave in the Material 1 I’eaches the gauge number 
1 . 0 ^ remains oonstant (with a eonstant stress magnitude) only 

during the first interval (or about 160 pseo) when the stress 
wave goes and returns between two ends of the stopper. Then 
the stress wave is transmitted and refloeted multiply at both 
ends of the stopper and starts to behave in a complex pattern. 
0 ^ is an almost step-like stress wave if the collision between 
a weight and a stopper is perfect, as both of them have homo- 
geneous cross sections normal to the longitudinal direction. 
However, as Is shown in Figures 6 and 7, does not rise 
sharply as step-like, indicating an imperfect collision. The 
pulse is actually a ramp type with a pulse rise time of about 


50 psec. Furthermore, as Is obvious in Figure 4, it is only 
until about 100 pseo counting from the wave front of when 


only (3^ is observed at the location of a gauge number 1. After 


this Interval, 


0 . 


will be superimposed with Its reflection at an 


predict the pattern of >r, ^ any longer. Even if ««, ^ is oonstant 
in time, the output from the No, 1 gauge will vary . Indeed as 
recorded in Fig. b and 7, the strain signal from the No. 1 gauge 
Increased after 100 ."see. This Increase must bo due to the super 
imposition of reflected waves from the adhesive Junction. From the 


rate of increase , one may infer that ‘U remains roughly constant 
until 160 sec. If the shear layer has sero thickness or its 
stiffness a is infinite, and the 


compound material Is assumed as a uniformly structured single 
body, variations of stress waves may be **»’eated by the theory of 
elasticity stress wave propagation In one diffienolon C63 . For 
our experimental specimen gf,j, (transmitted otress wave) * % 

0 ,^ (reflected stress wave) * % o^. Dot-dash linos In Figures b 
and 7 are for predictions obtained from the one~dlmenslonal wave 
propagation theory, assuming a uniform structure for the adhesive 
Junction, In this computation, the incoming stress wave is taken 
to be a ramp type pulse with a rise time of 30 usee. By comparing 
theoretical results in one dimensional wave theory with experi- 
mental results, one finds that results from the specimen A with 
a very thin adhesive layer agree well with theoretical predictions. 

Prom comparison of two types of specimens A and B, one recog- 
nises the significant difference in values of no.a-S. While the 
strain observed in an A *,’50cimen is consistent with that of one 
dlinenolonal theory, it hardly propagates in a B specimen that has 
a thick shear layer. Apparently, the diffoi’once is attributed to 
the thickness of shear layers. The magnitude of the shear stress 
transmitted to the Material 2 from the Material 1 is inversely 
proportional to the thickness b. For small b, the transmitted 
shear stress is large, and two materials separated by an adhesive 
layer form nearly uniform single body structure. Even though the 
A specimen hints the similar tendency, we found that the strain 
measured at No. 2-1 in the material, particularly in the B speci- 
men, is larger than not only the theoretical prediction, but the 
incoming stress wave itself. If the shear layer is thinner and 
the whole adhesive Junction can be regarded as a single body, then 
wave reflection takes place at the discontinuity where cross section 
changes and the wave magnitude passing through no. 2-1 have to be 
smaller than the incoming wave. However, experimentally the 
strain at No. 2-1 is larger than the incoming wave. For the thick 
shear layer the shear stress acting on the layer is small. It is 
not transmitted to the material immediately, even after it 
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pi»opa6atea through the Material 1 and reaches the adhesive 
junctions. While the incoming stress wave stays in the Material 
1 inside of the junetlon, it transmits the stress to Material 
V.hen the force i . transmitted to Material 2 via shear stress, 
Material 1 also receives its reaction and the gauge no, 3-1 
installed deeply in the junction records the strain (stress) 
larger than not only theoretical predicted transmitted wave, 
but the Incoming stress wave g^ itself. That is to say, the 
reflection plane for stress waves shifts effectively inward. If 
the shear layer gets thicker, the effective reflection plane moves 
further Inward. Accordingly, an incoming wave penetrates into 
Material 1 inside of the junction and then is superimposed with 
its reflections . Then It is within an adhesive junction where 


the first largo stress acts on Material 1. If Material 1 breaks 
down a little more than the Imeoraing stress, then the shear break- 
down of fibers could happen inside of this junction. At any rate, 

tnOV’i e* *t *1 mri** K rk m t » ^ »l ^ ■!» 4« jsw ^ ^ *-4 4 W 

vsivv/A ^ utr oUi i uO^:? 

stress far from the layer within the junction. But it fails to 
explain the details, much less the shear stress acting on a shear 
layer, that could be examined only by numerical analysis. 


5.2 Numerical analysis results 

As Is clear from Figures 6 and 7, the stress caused by the 
impulse of a fallen weight in a pull-out material (fiber) may be 
approximated in terms of a ramp type pulse with a pulse rise 
time of 50 psec. Therefore, we have examined in numerical anal- 
ysis two types of loadings, one with a ramp type pulse of 50 ysec, 
and another with an extremely sharp step pulse. We used the /SB 

following parameters ! 

2 

Young's modulus (for Al) E = 7000 (kg/mm ) 

Density (for Al) r = 2.75 (g/cm^) 

p 

Gravity acceleration g = 9.8 (m/s ) 

>-) 

Shear elasticity (for Araldlte) G = 50 (kg/mm*") 
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We assumed a lumped mass matrix of Equation (13) and a con- 

. . . . i » «r*. _»• , #«►__ J.* J» 


numerical method was employed to integrate equations of motion. 
Equation (8), while a unit interval for time integration is 
taken to be about 1/3 of the «ime required for the stress to 
pass through a finite element. Taking into consideration the 
phenomenon of wave propagation, we have divided the model in the 
figure into finite elements uniformly. We have used a coraputirr* 
HITAC 8300 (with memory capacity of 1^13 kb) . 


g . a . 1 


In order to make a valid comparison with experim?-nts, com- 
putation was done for the model in Figure 8a, which approximates 
the actual experimental specimen. Stress of a ramp type pulse 
with a rise time T ■ 30 ysec is assumed to act on the pull-o .t 

side of the model. The magnitune of the stress is taken to be 
P 

0.5 kg/mm . The length of an element is 2 mm for this computa- 
tion, We analysed the pull-out phenomenon as variations of 
strain at various sites in the model, until the reflected wave 
arrives there either from the loading end or from another free- 
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Pifture 8. Finite element model for calculation 
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Figure 9. Calculated results of strains and shear stress 
in shear layer for A specimen: In case of ramp type pulse 
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Figure 10. Calculated results of strains and shear stress 
in shear layer for B specimen: In case of ramp type pulse 
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Results are shown in Figures 9 and 10 for the shear layer 
thickness b » 0,05 mm and b « 1 mm, respectively. The time t 
is clocked from the instance the load is applied impulsively at 
the loading end. Prom the numerical curves agreement is remark- 
able qualitatively as well as quantitatively with experiments, 
except for the gauge No. 2-1 for the shear layer of b = 0.05mm. 

It could explain even smaller patterns at No. 2-2, that is, the 
effect of a shear deformation layer which could not be understood 
in one-dimensional stress wave propagation theory with afore- 
mentioned single body hypothesis. The largest shear stress 
exists at the pull-out end, regardless of the shear layer thick- 
ness. For a pulse of a long rise time 50 psec, the shear stress 
Increases at the arrival of stress wave and when the stress 
remains constant it also stays constant equal to the one under 
static loading. Small fluctuations are observed only towards 
the end of a ramp pulse where the Incoming stress wave vanishes 
discontinuously-. These oscillations get slo’wer for a thicker 
shear layer. Oscillations are more visible under a step pulge 
loading, as will be discussed below. 

5.2.2 Step pulse 

The most interesting phenomenon of the Impulsive pull -out is 
found under a step pulse loading where an incoming stress is 
applied more abruptly than the ramp type pulse in order to approx- 
imate a literally impulsive loading. Numerical analysis was 
made for the case where an incoming step pulse is realized with 
an extremely quick rise time T = 1 ysec on one end of the pull- 
out material (fiber). The model used in the analysis is shown 
in Figure 8b. To obtain an approximate solution for wave qua- 
tions Equations (1) and (2) for a sharply rising pulse, using 
the finite element method one has to solve an equivalent oscill- 
atory equation. Equation (1^0, where all the masses of waves 
with infinite number of degrees of freedom are concentrated on 
every nodal point that counts finite. For a better approxima- 
tion, for a sharper pulse one has to use as fine a division for 
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Figure 11. Calcualted results Figure 12. Calculated results 

of strains and shear stress in of strains and shear stress in 

shear layer for A specimen: in shear layer for B specimen: in 

case of step pulse. case of step pulse. 






eaah element an possitile. We solved the set of differential 
equations numerieally using the Runge-Kutta-Glll method. Prom 
consideration of the computer memory capacity and In order to 
minimise the number of elements to avoid a large accumulation 
of errors In the numerical method, we took the e3ement length 
of It " 1 mm. Hhus, one can examine only a shore time interval 
In the pull-out phenomenon, yet this Interval is long enough to 
understand the transitional phenomenon. Figure 11 shows numer- 
ical results for a b » 0.05 mm thick shear layer, One may notice 
a few over chutes at the top of the strain wave at the no. 1 cite 
sensitive to an Incoming stress wave. These may be unavoidable 
errors in the finite element method for wave equations. Although 
they cannot be ignored for a detailed study, they are negligible 
for an understanding of a global picture of the pull-out pheno- 
menon. The Inter-lamlnar shear stress fluctuates In time. The 
fluctuation is much less visible in response to a ramp type pulse. 
The fluctuation in inter-laminar shear stress finally decays to 
the value under static loading. On the other hand, the impulsive 
pull-out shear stress reaches a 20^ larger value over the static 
case in a transitional response period. Also, the stress varia- 
tion along a fiber is more complex. One nu.y assume that the fluc- 
tuation of the shear stress is generated in the adhesive Junction 
in the following mechanism. Firstly, an injected stress wave 
propagating in a fiber transmits the force towards the outer 
material, and pulls it abruptly through the shear layer, but the 
outer material cannot respond to it due to its inertia and then 
the shear stress increases. Within a fiber, particles move along 
the palled direction first, but then they are decelerated. Nextly 
the outer material starts being pulled, the relative displacement 
between the material and the fiber gets smaller and the shear 
stress also decreases . Particles in the outer material, by iner- 
tia, keep moving in the pulled direction, and the shear stress 
reduces further. Then the decrease of the stress forces the par- 
ticle deceleration that is caused by the strong shear stress 
decz’ease. Consequently, the shear stress increases again in the 
shear layer. This whole cycle repeats many times and there is a 
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fluctuation of the Inter-laminar shear stress. Therefore, the 
fluctuation must depend on the thickness of a shear layer. As 
is seen in Figure 12, it is less frequent for a 1 mm thick 
shear layer. This frequency is quite Insensitive to an element 
length as long as the variation of the division size for finite 
elements is sufficiently less frequent. This insensitivity could 
be characteristic of the Impulsive pull-out phenomenon. 

6. Conclusion 


Numerical analysis and experimental results are presented on 
the pull-out phenomenon in fiber-enforced composite materials at 
the high rate of loading accompanying impulsive breakdowns. The 
finite element method was used assuming an inter-laminar shear 
layer as the interface between a fiber and the matrix. If the 
fiber is pulled out Impulsively, the shear stress is time depen- 
dent and decays into the value under static loading. Under step 
pulse loading, the inter-laminar shear stress fluctuates, finally 
decaying into the static value. The maximum shear stress under 
the impulsive loading assumes a value more than 20^ larger than 
the static value . 


(A lecture given at the 5th Symposium on PRP by the Japan 
Society of Material Science, May 31, 1979). 
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